Central peak position in magnetization loops of high-T c superconductors 
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Exact analytical results are obtained for the magnetization of a superconducting thin strip with a 
general behavior J C (B) of the critical current density. We show that within the critical-state model 
the magnetization as function of applied field, B a , has an extremum located exactly at B a — 0. 
This result is in excellent agreement with presented experimental data for a YBa2Cu307_^ thin 
film. After introducing granularity by patterning the film, the central peak becomes shifted to 
positive fields, B cp > 0, on the descending field branch of the loop. Our results show that a positive 
B cp is a definite signature of granularity in superconductors. 



PACS numbers: 74.25.Ha, 74.76.Bz, 74.80.Bj 

The investigation of magnetic hysteresis loops (MHLs) 
is a widely used tool to characterize superconducting 
samples, and in particular, to estimate the critical cur- 
rent density and its dependence on magnetic field. An 
ever-present feature of the MHLs is a peak in the mag- 
netization located at an applied field B cp near zero. The 
central peak is formed due to a field dependence of the 
critical current density, J C {B), monotonously decreasing 
at small fields. When the sample is a long cylindrical 
body placed in a parallel applied field, the peak position 
can be calculated analytically within the critical-state 
model for several J C (B) dependences One always 
finds on the descending field branch that B cp < 0, and 
similarly B cp > on the ascending branch of large-field 
loops. In simple terms the shift is a consequence of the 
local flux density, B, lagging behind the applied field. 

Also experimentally there are abundant observations 
of negative B cp on the descending field branch. How- 
ever, in some cases one finds B cp very close to zero or 
even shifted to the positive side so that the peak occurs 
before the remanent state is reached |^|6| ■ The explana- 
tions for this behavior are controversial. On one hand, 
a decreasing sample thickness is known from experiment 
to shift the peak toward B cp — [0J|]. This is also in 
agreement with recent numerical results obtained for the 
critical-state model with B-dependent J c |^,^|. However, 
the numerical calculations have not predicted sufficiently 
large shifts to bring the peak to the anomalous side of 
the loop. 

A different explanation of the shift in B cp was sug- 
gested for granular materials, jlO| where a demagneti- 
zation effect of the grains can be important. The mag- 
netization of a granular superconductor is determined by 
both intra- and inter-grain currents. The latter represent 
large-size current loops, and can give a major contribu- 
tion to the magnetic moment. These inter-grain currents 
are essentially determined by the magnetic field just at 
the grain boundaries. Those local fields are, in turn, 
strongly influenced by the intra-grain currents, and can 



vary ahead of the local B expected without granularity. 
As a result, the central peak is encountered earlier in 
the MHL. Models allowing quantitative MHL calcula- 
tions with account of the grain-induced demagnetization 
effect [^|JTl]l successfully reproduced the peak at B cp > 0. 

Although the significant influence of both sample 
shape and sample granularity on the MHL is generally 
recognized, the interplay between them has never been 
addressed. In particular, it is not known if a positive B cp 
is accessible for thin enough samples, or if B cp > on the 
descending field branch is a definite signature of sample 
granularity. The present work aims to answer this ques- 
tion. We show first analytically that for any J C (B) the 
central peak is located exactly at B cp — for a thin uni- 
form strip, a result we also confirm by experiment. The 
key role of granularity is then demonstrated by measure- 
ments on an artificially granular thin film, where we find 
B cp > 0. 




FIG. 1. Superconducting strip in applied magnetic field. 

Consider a long thin superconducting strip with edges 
located at x = ±10, the y-axis pointing along the strip, 
and the z-axis normal to the strip plane, see Fig. [I]. 
The magnetic field, B a , is applied along the z-axis, so 
screening currents flow in the y-direction. Through- 
out the paper B is the z-componcnt of magnetic induc- 
tion in the strip plane. The sheet current is defined as 
J(x) = J j(x, z) dz, where j(x,z) is the current density 
and the integration is performed over the strip thickness, 
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d < w. 

From the Biot-Savart law for the strip geometry, the 
flux density is given by pa 



B(x)-B a = -^ J 



Mo f w J( u ) du 



(1) 



Assume that the strip is in a fully penetrated state, i.e., 
the current density is everywhere equal to the critical one, 
J(x) = sign(x) J c [B(x)]. This state can be reached after 
applying a very large field, and then reducing it to some 
much smaller value. The field distribution then satisfies 
the following integral equation, 



B(x)-B a = -«> 

T Jo 



J c [B(u) 



udu . 



(2) 



In the remanent state, B a = 0, the flux density pro- 
file B(x) has an interesting symmetry. This is seen by 
changing the integration variable in Eq. (Q) from u to 
v = V 'w 2 — u 2 . We then obtain 



B{x) — B a — — / - - - — vdv. 



(3) 



Substituting x — > \Jw 2 — x 2 into Eq. (Eh one obtains 



similar equation for B{\/w 2 
B(Vw 2 -x 2 )-B a = -^ \ 



Jc[B{u)] 

■t. 7, ~ udu . (4) 



w — x A — u* 



By comparing equations (||) and (]!]) at B a = we con- 
clude that 



B(x) = -B(^/w 2 -x 2 ) 



is generally valid if J c depends only on the absolute value 
of the magnetic induction, J C (|_B|). This symmetry is im- 
mediately evident for the case J c = const., i.e., the Bean 
model, when Eq. (0) reduces to 



B(x) = B a + B c In 



Vw 2 — x 2 



HqJc 



(6) 



In a similar way one can prove also another symmetry 
relation valid at B a — 0, namely 



V(x) = V(\/w 2 - x 2 ) , 



(7) 



where V(x) = dB(x)/dB a . 

Consider now the magnetic moment per unit length of 
the strip, M = 2 L J(x)xdx. Differentiating M with 
respect to B a and taking into account that B(x) changes 
sign at x — w* — w/y/2 one has after splitting the inte- 
gral into two parts 



dM 



dJ c (\B(x)\) 




V(x)xdx. (8) 



Then, replacing in the second integral x by y/w 2 — x 2 
and using Eq. dq), we come to 



dM 

dB~a 



at B a = 



(9) 



Consequently, on a major MHL the magnetic moment 
has an extremum in the remanent state for any J C (B)- 
dependence. This is illustrated in Fig. || showing the 
central part of MHLs obtained by an iterative numeri- 
cal solution of Eq. (||) for three J C (B) dependences. The 
M(B a ) extremum depends on the type of J C (B) depen- 
dence. It is a maximum for a decreasing J C (B), and a 
minimum if J C (B) has a pronounced second peak, the 
so-called fishtail behavior. 
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FIG. 2. Central part of descending branches of the MHL 
for different J C (B) dependences, J c /J c o = F(\B\/B c o)- (1) - 
T{b) = (1 + 2.56 2 )e- 62 / 4 ; (2) - T{b) = (1 + 0.256)- 1 ; (3) - 



(5) T(b) = e~ b/2 . Here B c0 = u J c0 /ir, Mo = w 2 J c0 . 



Two YBa2Cu307_a epitaxial thin films were chosen for 
a comparative study; (A) a uniform film of regular shape, 
and (B) a film with artificial granularity. Sample (A) 
was prepared by laser ablation on an MgO substrate [^3) . 
The thickness was 200 nm, and the sample was patterned 
by chemical etching into a rectangular shape of dimen- 
sions 0.66 x 1.4 mm 2 . The homogeneity of the sample 
was tested by magneto-optical imaging, where it showed 
complete absence of any visible defects JlJ] . 

Sample (B) was laser ablated to a 150 nm thickness on 
a LaA103 substrate. The film was patterned by means of 
electron beam lithography into a hexagonal close-packed 
lattice of disks with 2r = 50 /im diameter. The disks are 
touching each other at the circumferences in order to en- 
able the flow of inter-disk currents. The width of the con- 
tact region is e = 3.5 /im, see Fig. ||. The superconduct- 
ing disks can be considered as grains and the connections 
between them as inter-grain junctions. The overall size 
of the sample was 4x4 mm 2 , comprising ps 8000 disks. 
The transition temperature after the patterning process 
is 83 K. 
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FIG. 3. Polarization image of the sample (B) with artificial 
granularity. The black disks are superconducting and have a 
diameter of 2r = 50 /im. The width of the contact region 
between the disks is e = 3.5 /xm. 



ferent J C (B) dependences. There is evidently here a pro- 
nounced central peak at all temperatures. Within the ex- 
perimental resolution the position of the peak is located 
at zero applied field. The observation that the peak re- 
mains at B a = over the entire temperature range is in 
full agreement with our general analytical result. 

Note that the result was derived for an infinite strip 
fully penetrated by the magnetic field. In a thin sam- 
ple, in contrast to long cylinders, the flux fronts move 
in response to a changing applied field at an exponential 
rate |l2[ . Therefore, a fully remagnetized state is reached 
quickly after reversing the direction of a field sweep, and 
our experimental conditions are consistent with the as- 
sumptions made in the theory. The only exception is the 
finite length of the strip. However, our results show that 
it is not a crucial factor for the peak position. 



The magnetization measurements were performed us- 
ing a vibrating sample magnetometer (VSM) PAR Model 
155 with computerized control and data processing. The 
magnetic field is generated by a conventional magnet 
with -B a ,max ±2 T. The field was applied parallel to the 
c axis, which is perpendicular to the film plane. The 
magnetization loops were always measured after zero- 
field cooling of the sample. The maximum applied field 
for all measurements was significantly higher than the 
field giving a fully penetrated state, as verified by the 
magneto-optical method. To focus on our main issue, 
only the data measured on the descending field branch 
of the MHLs near B a = are presented. 
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FIG. 4. Central part of descending branches of the MHL 
for a thin YBCO strip for different temperatures. The peak 
in magnetization is located always at zero applied field. 

Figure ^ shows the central part of the MHLs for the 
uniform film (A). The measurements were carried out 
over a wide range of temperatures in order to probe dif- 



0.04 



E 

CD 



0.02 



0.00 



5 K 



15 K 



"□□□□□ Dm 20 K 



30 K 



40 K 
50 K 



0.0 



0.1 0.2 
B a (T) 



0.3 



FIG. 5. Central part of the descending branch of MHLs for 
the sample (B) with artifical granularity. The solid line indi- 
cates position of the magnetization peak which is located at 
positive fields at all the temperatures. 

The MHLs for the sample (B) with artificial granular- 
ity are shown in Fig. |5|. The central peak displays here 
different behaviour as compared to the uniform film (A) . 
For all temperatures B cp > on the descending field 
branch. The shift of the central peak becomes more and 
more pronounced with decreasing temperature reaching 



B c 



+40 mT at T = 5 K. 



In this sample, the "grains" and their interconnections 
are of the same material, and hence have equal critical 
current density. Due to the patterning, the spatially av- 
eraged inter-grain current density J' r is less than the cur- 
rent density, J-T , of the superconducting material itself. 
For the hexagonal pattern, see Fig. [| one has J* r = r]J-! r , 
where r\ depends on the direction of J' r , within the limits 
e/2r < r\ < e/^/ir. Futhermore, the magnetic moment 
M 9r due to intra-grain currents and moment M tr due to 
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inter-grain currents are related as 

M ar /M tr = (Jf /Jl r ){r/R) « 2r 2 /(eR), 

where 7? is the overall size of the sample. In the present 
case this ratio is 1/6, showing that the inter-grain contri- 
bution to the total magnetic moment is dominant. A con- 
siderable grain-induced demagnetization effect is there- 
fore expected for this sample. Its granularity, being the 
only essential difference from the uniform film, is the only 
possible origin of the observed shift to B cp > 0. 

The superconducting Bi-2223 tapes are known to have 
a granular microstructure and the two contributions to 
the magnetic moment are of comparable magnitude |Q. 
The grain demagnetization effect has been suggested as 
an explanation for the observed shift of the peak to 
B cp > on the descending MHL branch in these ma- 
terials Our present results give strong support to 
this physical picture. 
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FIG. 6. A drawing illustrating the main conclusion of the 
paper. Long samples in parallel magnetic field show the peak 
on the descending field branch of MHL at a negative field. As 
the sample height decreases the peak is shifted towards zero. 
It is exact zero in the limiting case of a uniform thin strip. A 
granular thin strip has the peak at positive fields. 

In conclusion, we arrive at a following general scenario 
concerning the central peak position, see Fig. ^. Long 
samples in parallel magnetic field show in their MHLs a 
central peak at a negative applied field, i.e. after pass- 
ing through the remanent state on the descending field 
branch. As the sample thickness decreases the peak po- 
sition is shifted towards zero. It is located exactly at 
^cp = 0, in the limiting case of a uniform strip of in- 
finitessimal thickness. Granularity always leads to a shift 
of -B cp in the positive direction on the descending field 
branch. Thus, for a granular thin strip the peak is lo- 
cated at a positive field, i.e., before the remanent state is 
reached. The origin of this effect is that granularity in- 
duces demagnetization fields which strongly modify the 
inter-granular currents via their independence. 
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